Magnetization Process of Single Molecule Magnets at Low Temperatures 
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Magnetic clusters, such as Fes and Mni2, that make 
up the core of large organometallic molecules, behave 
at low temperatures as large single spins. Accordingly, 
these molecules have come to be known as single molecule 
magnets (SMM's) 0. In crystals, magnetic anisotropy 
energies U inhibit magnetic relaxation of SMM's, which 
can consequently take place at very small temperatures 
only through magnetic quantum tunneling (MQT). Dipo- 
lar interactions play then an essential role. They can 
give rise, upon tunneling, to Zeeman energy changes of 
nearly 1 K. This exceeds by many orders of magnitude 
the ground state tunnel splitting energy A that would 
follow from perturbations by higher anisotropies for Feg 
and Mni2 |3|- Energy conservation would make MQT, 
which has been observed experimentally impossible 
for the vast majority of spins in the system. Hyperfine 
interactions between the tunneling electronic spins of in- 
terest and nuclear spins open up a fairly large tunneling 
window of energy Ew such that tunneling can occur if the 
Zeeman energy change 2eh upon tunneling is not much 
larger than 0. More precisely, the tunneling rate F' 
for spins at very low temperature, is given by 



(1) 



where F is some rate (whose value is not important for 
our purposes), ri{x) ~ 1 for | x |< 1, ri{x) ~ for a; > 1, 
and Eiu > A. Other theories for MQT of SMM at very 
low temperatures have also been proposed 'ij . We adopt 
Eq. here, regardeless of theory or physical mechanism 
behind it. We let ri{x) = 1 for | a; |< 1 and ri{x) = for 
a; > 1, and refer to £„, as the tunnel energy window. 

The interesting early time relaxation 1 — m cx y/t 
of an initially magnetized system has been predicted 
0, observed experimentally @, further explained [tII, 
and widely discussed 0. An unpredicted related phe- 
nomenon was later observed by Wernsdorfer et al. 0: 
the magnetization m of a system of Feg SMM's increases 
as i/t, where t is the time after a weak magnetic field 
is applied to an initially unpolarized system. Conve- 
niently, this latter effect seemed to be independent of 



system shape. Interesting questions arise: is this a uni- 
versal effect to be found in all MQT experiments? If not, 
what does it depend on? How many time decades does 
the y/t regime cover? What is the final steady state mag- 
netization? No explanation or simulation that we know 
of has been offered. We address these questions here. 

We report Monte Carlo (MC) results that reproduce 
the m (X Vt behavior of initially unpolarized systems. 
We show that this arises from correlations that develop 
between spins and local magnetic dipolar fields, while 
cooling to low temperatures, before finally quenching to 
experiment. Furthermore m{t) depends on the cooling 
protocol only through the final energy — reached just 
before quenching. 

The main results obtained follow. All energies and 
magnetic fields are given in terms of and hd, respec- 
tively, where —2ed is the energy of two S spins that lie 
on sites a distance away, that point along the line join- 
ing their two sites, a is the side of a cubic unit cell, 
hd = £d/ {g^J'BS), g is the gyromagnetic ratio, and 
is the Bohr magneton. We also let a stand for the rms 
value of h for a disordered spin configuration • After 
quenching and applying a field < 1 at t = 0, 



m{t) ~ bSaEwHa ^F(Fi, (t/Eu), o'/Ziq), 



where h • 
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for Tt < 1 

for 1 <rt < (cr/e^)i/P 
for (a/e^)i/P < Ft, 
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where, ho = 2(27r) /3 ' pu, py is the number of spin 
sites per unit volume, p ~ 0.5, for simple cubic (SC) 
lattices, and p ~ 0.7 for body centered cubic (BCC) 
and face centered cubic (FCC) lattices. These results 
are inferred from MC simulations in which the energy of 
the magnetic system is held constant in time as well as 
from arguments given below. For magnetic systems that 
readily exchange energy with the lattice, for which the 
energy is not a constant of time, results are briefiy dis- 
cussed in the closing remarks. We first describe the sim- 
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FIG. 1: (a) m/e-w versus Vt for an applied external field H — 1 
for the shown values of £„. First, thermalization took place 
from an initially disordered configuration up to the time when 
Ea reached the value 0.58. The dashed, continuous, and dot- 
dashed lines are for curves calculated from Eq. IIH , for — 
0.05, 0.2 and 1, respectively, (b) Same as in (a) but for the 
shown values of Sa and Sw = 0.1. The continuous line, for 
the solution from Eq. (Illll . (c) Same as in (a) and (b), but 
for different lattices, with = 0.05. For SC, BCC, and FCC 
lattices, previous partial thermalization took place at Ta — 10, 
20 and 60, till Sa = 0.31, 0.28, and 0.53, respectively. The 
continuous, long dashed, and short dashed lines are for curves 
calculated from Eq. itTTl for FCC, BCC, and SC lattices, 
respectively. There are no adjustable parameters. 



ulations. We use the MC method to simulate magnetic 
relaxation of Ising systems of ±S spins, on simple cubic 
lattices with periodic boundary conditions, that inter- 
act through magnetic dipolar fields and flip under rules 
to be specified below [ll|. The system is first allowed 
to evolve towards thermal equilibrium at some "high" 
temperature Ta- We assume ksTa ^ U/10, which im- 
plies that spin reversals then take place mostly through 
classical thermal processes. Accordingly, spin flips are 
then governed by detailed balance rules, and Eq. ^ is 
not enforced. For reasons that are given below, we also 
impose the restriction Ta > Tq, where Tq is the long- 
range ordering temperature. One may think of this first 



process as a waiting stage that the systems may have 
to undergo in the cooling process, before quenching to 
a lower temperature where a tunneling experiment (as 
in Ref. ^ can be later performed. Let this first stage 
end at t — with sudden cooling of the system to a 
"very low" temperature, that is, to a temperature be- 
low, roughly, 0.2U/{SkB) Hill- Accordingly, Eq. Q is 
then enforced on all spin flips for t > 0. As for detailed 
balance, we then proceed as follows. We assume that 
thermalization of a SMM system with the lattice does 
not take place (i.e., the energy is constant) at very low 
temperatures (but see Ref. [l3)- We meet this condition 
by enforcing detailed balance but using an appropriately 
chosen pseudotemperature Tu- [From an expression be- 
low Eq. ®, ksTu ~ (T^/2ea. Note that T^ > Ta, since 
— Eq cannot be smaller than the equilibrium energy at 
Ta]. We have checked that the mean energy is indeed 
constant under this rule. We do not report here results 
we have obtained applying detailed balance rules with 
T < Ta (applicable to systems where thermal relaxation 
to the lattice takes place ^3); but make a comment on 
them in the closing remarks. MC results for the time evo- 
lution of mjew in SC lattice systems, after a fleld H = 1 
is applied upon quenching, are shown in Fig. 1(a) for 
various values of e^,. Before quenching, the system was 
thermalized at Ta = 10 /ks for some time till the energy 
per spin reached the value —0.58. Clearly, m scales with 
up to a crossover time of, roughly, lOF^^/e^, where 
TO levels off. Within the time range 1 < Ft < 10/e^, 
TO cx ^/t. Monte Carlo results that show how m scales 
with Ea are exhibited in Fig. 1(b) for = 0.1. Note 
also that TOe, the leveling off value of to, scales with e^, 
and, as argued below, it scales with ct"^ as well, in agree- 
ment with Fig. 1(c). Results for different cubic lattices 
are shown in Fig. 1 (c) . The log to versus log t slopes in 
the intermediate time regime are clearly lattice structure 
dependent. Much of Eqs. (2-5) is inferred from these 
graphs. 

For most of the rest of the paper, we try to understand 
these results. Let's first examine the physics of the wait- 
ing stage. We assume the system is first cooled to some 
temperature Ta that is above the ordering temperature 
Tq, but not infinite. We also assume that /cbTq > U/10. 
It then follows from Arrehnius' law, t = tq exp{U/kBT), 
that over barrier spin flipping readily takes place in the 
laboratory within a second's time if tq ^ 10~^ s. Some 
correlation between spin Si and field hi at each site i 
can therefore be established, but no long-range order can 
develop if Ta > Tq. Assume that either Ta 3> To or 
that the time spent in the waiting stage is so short that 
the probability density function (PDF) p{h) that a ran- 
domly chosen site have field h is reasonably approximated 
by p{h) (X exp(— ft,^/2cr^). On the other hand, the con- 
ditional PDF to find ztS given a field h acting on the 
spin fulfills, in equilibrium, p(±S' | h) cx eKp(±h/kBT). 
Now, since the joint probability density p{±S, h) that. 



on a randomly chosen site, one find h acting on ±S is in 
general given by p{±S, h) = p{±S \ h)p{h), 

p(±5,/i)oce-(''T^'/'=-^)'/2-', (6) 

follows in equilibrium. Therefore, the mean energy is 
—a^/2kBT. The replacement cP' jlkBT ^ Sa-, generalizes 
the above equation to 

for all times up to equilibration. Then, to first order in 

PT {h) - vv {h) ~ teafT-3e-''V2-^ ^ (8) 

where P'\{h) — p{S, h) and P[{h) = p{—S, h). All points 
for — pi{h)]/{hea) obtained from MC calculations 

collapse onto a single curve in Fig. 2a, in agreement 
with Eq. ||SJ|. Deviations, of higher order in Ea/c, from 
Eq. © do occur. They are within 10% even if complete 
thermal equilibration is allowed to take place as long as 
T > 10, i.e., above approximately ATq. 

We now examine the time evolution of the system af- 
ter abruptly cooHng it, at time t = 0, to a tempera- 
ture below, roughly, Q.2U /{Sks)- Then, spin flips up to 
I |< S states can be neglected, and tunneling through 
the ground state doublet is the only available path for 
spin reversals. Accordingly, spin flips are allowed only if 
the spin's Zeeman energy is within the tunnel window. 
A field H is applied for all t > 0. Now, if either the 
system is in thermal contact with a reservoir at a tem- 
perature such that /cfiT ^ e^, or the energy is constant 
and sufficiently high such that ksT^ 3> ew^ then 

m :^2r J dh r]{H + h)f{h, t), (9) 

where f{h, t) = pi{h, t) - p^{h, t). 

It is worth pointing out that — /(ft., 0) is given by Eq. 
©, and therefore /(ft,0) = -{2hea/ (j'^)p{h,Q) at t = 0, 
where p{h, 0) is the PDF regardless of spin orientation, 
i.e., p{h) = pi{h) However, f{h,t) oc hp{h,t) 

does not hold for t > 0. This point is illustrated in Fig. 
2(b), where MC results for both f{h,t) and p{h,t) are 
shown for some non-zero times. _p{h, t) is invariant for 
times ^ r. As reported in Ref.[lJ|, a sharp dip, of ap- 
proximately e^, half-width, does develop in p{h), but only 
at much later times, and then not when ksT ^ e^. On 
the other hand, a hole does show up in f{h,t) in Fig. 
2(b), as in Wernsdorder's experiments 0, performed at 
T = 40 mK (which is roughly 10 times as large as Sw/kB) 

The time development of the hole in f{h, t) is illus- 
trated in Fig. 2(c). Note that the hole deepens, but 
its width remains approximately constant while Ft ^ 1, 
and then, under the conditions stated above Eq. lO, 
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FIG. 2: [p^{h)~pi{h)]/{hea) for a system of 16 x 16 x 16 spins 
at Ta = 10 (recall To ~ 2.5), before quenching. The shown 
temperatures are given in units of £d/fcfl. The system was 
cooled to these "high" temperatures from an infinite tempera- 
ture, and allowed to thermalize for a time, till the shown ener- 
gies were reached. AH points stand for averages over 1.5 x 10^ 
runs, (b) [pi(h) - p^{h)]/ (hea) and p(h) = p^(h) +PT{h) at 
times t = 0, Ft — 0.62 (discontinuous line), and Ft — 2.19 
(continuous line), for the same system as in (a), after the sys- 
tem was first thermalized at T = 10 till Sa — 0.25, was then 
further cooled to T = 1, and H = 2 was then applied. A 
tunnel window Ew =0.1 was enforced. All points in (b) stand 
for averages over 4 x 10* histories, (c) Same as in (b) but for 
Ft = 0.5, 1, and 4; in addition, while t < 0, the system had 
been partially thermalized at T = 10 till Sa ~ 0.247. 



/ = — 2r/. Therefore, since m equals the area coverd by 
the hole [i.e., m{t) ^ J dh [/(ft,0) - f{h,t)]], 

m~2£^/(-i/,0)(l-e-2rt) (lo) 

if Ft < 1 and < a. Using Eq. ©, Eq. (3) follows. 
The value rUe that m{t) levels off to after a sufficiently 
long time, that is, Eq. (5), can be estimated as follows. 
For Ft > 1, f{-H,t) < /(-i/,0), and the hole only 
becomes broader, but it cannot become wider than the 
field distribution p{h). The final area covered by the hole 
is therefore ~ 2f{—H,0)a, which is the estimated value 
of TOg, in rough agreement with the expression for rUe 
below Eq. (5). 
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More detailed considerations underlie Eq. (4). We 
have derived the equation, 

a;~ciW--C2/ dT— , , (11) 

V TT Jo ^{t -T) + 

where, ci and C2 are constants to be specified below, 
uj{t — t) is the inverse of the PDF that the field h on a 
randomly chosen site be the same at times t and t. Quan- 
tity m follows from Eq. Hll() letting m — x£a£wH / , 
ci — 4, and C2 = 2. uj{t — t) depends on i — r through 
the probability <f>[t — t) that a spin point in opposite di- 
rections at times t and t. To make progress, we make 
the approximation oj ~ rnh\[2'Kho4>,2a^/T:^/(p\. The ap- 
proximation for (/) <C 1 follows from the Lorentzian PDF, 
of half width /lo^, that ensues when a small fraction 
of sites are randomly occupied by spins [13. (The two 
factors of 2 come from the fact that flipping an already 
present spin S is like placing a 25* spin on an unoccu- 
pied site.) The approximation for 0^1/2 follows from 
the Gaussian distribution that holds then. In between, 
the interpolation checks, within some 10%, with our MC 
results. Finally, 0(t — r) must be evaluated. For this 
purpose, an equation for the fraction of spins n{t) that 
have flipped at least once before time t is derived [isj . 
It is Eq. (|ll|l using n = xe^/a, ci = 1 and C2 = 1. 
We then use (/) « n/2. The functional dependence of F 
shown in Eq. (2) follows by careful inspection of these 
equations. Numerical calculations yield the curves shown 
in Figs. 1(a), (b), and (c). The exponent p depends on 
lattice structure, through (t//io in the expression for a;(0) 
above. The agreement with our MC results exhibited in 
Fig. 1(c) is reassuring. 

A couple of final remarks follow. Equations (2-5) are 
for magnetic systems that do not exchange energy with 
a heat bath (i.e., the lattice, usually) at very low tem- 
peratures. Then, m{t) — > as t — > oo. When heat 
exchange does takes place readily, as in some systems in 
Ref. |l j . then our simulations show that m(t) eventu- 
ally crosses over from the value given by Eqs. (2-5) to 
mth{t) ~ 0.3{ew/a)^HTt, for H <1. This happens when 
iTith{t) becomes the larger of the two. Later on, nith lev- 
els off to a quasi steady state value (not the final thermal 
equilibrium value), which depends on both system shape 
and lattice structure, as expected. (If a field is not ap- 
plied inmediately upon quenching, but later, the quasi 
steady state value of mth is also affected.) The {ew/<j)^ 
dependence suggests that (see Ref. ^3); contrast with 
constant energy magnetization processes, magnetic or- 
dering takes place while the system magnetizes. 

Summing up, we have given MC and theoretically 
based evidence to show that the m (x y/t behavior ob- 
served in experiments on Feg clusters after quenching 
and applying a small field at t = is driven by cor- 
relations which are previously established in the system 
while cooling to very low temperature. We have estab- 
lished that the -\/t behavior is not universal. More gener- 



ally, m{t) (X tP, and p depends on lattice structure. The 
time range over which this behavior prevails, the value 
TOg that the magnetization later levels off to, and the 
crossover time to a final thermally-driven regime have 
been determined. More specifically, Eq. (2-5) have been 
inferred from MC simulations and Eq. Hll|) , and much of 
the relevant physics has been explained. 
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